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I. INTRODUCTION 

Many scientists have underlined the importance of a 
quantitative approach in social sciences [1-8]. In fact, 
statistical mechanics and social sciences have been al- 
ways linked to each other in a constructive way, due to 
the statistical character of the objects of study [9, 10]. 
On one hand, various discoveries, first made in the field 
of social sciences, introduced new concepts which turned 
out to be relevant for the development of statistical me- 
chanics and later of the science of complex systems. For 
instance, fat tails were found by Pareto in the distri- 
bution of wealth [11, 12]; the first description of finan- 
cial time series through statistical mechanics, made by L. 
Bachelier in his PhD thesis [13-15], also represents the 
first formalization of a stochastic process in terms of the 
random walk model; large fiuctuations were observed by 
Mandelbrot in the time series of cotton price [16]. On the 

\ other hand, physics has often represented a prototype for 
modelling economic systems. For example, many works 
of Paul Samuelson were inspired by thermodynamics; the 
analogies between physics and economics were studied 
by Jan Tinbergen in his PhD thesis entitled "Minimum 

'. Problems in Physics and Economics". Recent develop- 
ments of economics rely more and more on the theory 
of stochastic processes and the science of complex sys- 
tems [17]. 

The present paper considers some models of wealth 
exchange between individuals or economical entities, in- 
troduced independently in different fields such as social 
sciences, economics, and physics. We refer to them as ki- 
netic wealth- exchange models (KWEM), since they pro- 
vide a description of wealth flow in terms of stochastic 



wealth exchange between agents, resembling the energy 
transfer between the molecules of a fluid [2, 18, 19]. In or- 
der to maintain the discussion at a fundamental level, we 
limit ourselves to the following simple KWEMs: those 
introduced by Angle (A-models) [20-23], Bennati (B- 
model) [24-26], Chakraborti and Chakrabarti (C-mo- 
del) [27], and by Dragulescu and Yakovenko (D-mo- 
del) [28] . The goal of the paper is to discuss their general 
common features, formulation, and stationary solutions 
for the wealth distribution. We consider a heterogeneous 
KWEM, in order to illustrate how a simple KWEM can 
generate realistic wealth distributions. We also clarify 
some relevant issues, recently discussed in the literature, 
concerning the relaxation to equilibrium and the appear- 
ance of a power law tail of the equilibrium distribution 
in heterogeneous models. 

A noteworthy difficulty in the study of wealth or money 
exchanges based on a kinetic approach had been pointed 
out by Mandelbrot [29] : 

... there is a great temptation to consider 
the exchanges of money which occur in eco- 
nomic interaction as analogous to the ex- 
changes of energy which occur in physical 
shocks between gas molecules... Unfortu- 
nately the Pareto distribution decreases much 
more slowly than any of the usual laws of 
physics... 

The problem referred to in this quotation is that the 
asymptotic shape of the energy distributions of gases pre- 
dicted by statistical mechanics usually have the Gibbs 
form or a form with an exponential tail. The real 
wealth distributions, instead, exhibit a Pareto power law 
taU [11, 12, 30-32], 

fix) ^ llx'+'^ , (1) 
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with 1 < a < 2. However, it has become clear that 
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(a) the actual shapes of wealth distribution at interme- 
diate values of wealth arc well fitted by a F- or an ex- 
ponential distribution [21, 33 35], so that they can be 
reproduced also by simple KWEMs with homogeneous 
agents (see Sec. Ill); (b) KWEMs with suitably diversi- 
fied agents can generate also the power law tail of the 
wealth distribution [2, 19, 36] (see Sec. IV). This has 
opened the way to a simple, quantitative approach in 
modelling real wealth distributions as arising from wealth 
exchanges among economical units. 

The paper is structured as follows: In Sec. II a general 
description of a KWEM is given. In Sec. Ill the homo- 
geneous A-, B-, C, and D-models are discussed. Explicit 
analytical fitting forms for the equilibrium wealth distri- 
butions are given. In Sec. IV we discuss the influence 
of heterogeneity, taking the heterogeneous C-model as a 
representative example. In this respect, we analyze the 
mechanism leading to a robust power law tail. Some is- 
sues concerning the convergence time scale of the model 
and the related finite cut-off of the power law are dis- 
cussed. In Sec. V a unified reformulation of the A-, C-, 
and D-models is suggested, which in turn naturally lends 
itself to further generalizations. An example of gener- 
alized model is worked out in detail. Conclusions are 
drawn in Sec. VI. 



II. GENERAL STRUCTURE 

In the models under consideration the system is as- 
sumed to be made up of N agents with wealths {xi > 0} 
(i = 1,2,..., N). At every iteration an agent j exchanges 
a quantity Ax with another agent k chosen randomly. 
The total wealth X = Xi is constant as well as the 
average wealth (x) — X/N. After the exchange the new 
values x'j and x^ are {x'j , > 0) 

x'k = Xk + Ax . (2) 

Here, without loss of generality, the minus (plus) sign has 
been chosen in the equation for the agent j (k). The form 
of the function Ax — Ax{xj,Xk) defines the underlying 
dynamics of the model. 

In KWEMs, agents can be characterized by an ex- 
change parameter lo E (0, 1] which defines the maximum 
fraction of the wealth x that enters the exchange pro- 
cess. Equivalently, one can introduce the saving param- 
eter A = 1 — u), with value in the interval [0,1), rep- 
resenting the minimum fraction of x preserved during 
the exchange. The parameter cj (A) also determines the 
time scale of the relaxation process as well as the mean 
value (x) at equilibrium [37] . If the value of to (A) is the 
same for all the agents, the model is referred to as ho- 
mogeneous (see Sec. III). If the agents assume different 
values Wj (Aj) then the model is called heterogeneous (see 
Sec. IV). Homogeneous models can reproduce the shape 
of the F-distribution observed in real data at small and 



intermediate values of the wealth. For w < 1 (A > 0), 
they have the self-organizing property to converge to- 
ward a stable state with a wealth distribution which has 
a non-zero median, differently from a purely exponential 
distribution. Models with suitably diversified agents can 
reproduce also the power law tail (1) found in real wealth 
distributions. 

In actual economic systems the total wealth is not con- 
served and a more faithful description should be used. It 
is therefore interesting to observe how the closed econ- 
omy models considered here, in which J2i is constant, 
provide realistic shapes of wealth distributions. This sug- 
gests that the main factor determining the wealth distri- 
bution is the wealth exchange. 

When the variation of wealths is not due to an actual 
exchange between the two agents but the quantity Ax 
is entirely lost by one agent and gained by the other 
one, the model is called unidirectional. Furthermore, it is 
possible to conceive multi-agent interaction models, not 
considered here, in which a number M > 2 of agents enter 
each trade. Then the evolution law has the more general 
form X- = Xi -|- Axi, with i = 1,. .. , M, J2f^i Axj = 0, 
and the Axj depending somehow on the wealths x, of the 
M interacting agents. 

III. HOMOGENEOUS MODELS 

A. Al-model 

Here we consider the model introduced by John An- 
gle in 1983 in Refs. [20, 21], referred to as Al-modcl (a 
different model of Angle, the One-Parameter Inequality 
Process, referrcid to as A2-model, is consider in Sec. Ill B 
below). The A-models are inspired by the surplus theory 
of social stratification and describe how a non-uniform 
wealth distribution arises from wealth exchanges between 
individuals. 

The Al-model is unidirectional and its dynamics is 
highly nonlinear. The dynamical evolution is determined 
by Eqs. (2) with Ax given as 

Ax = eujiVjM Xj - (1 - T]j^k) Xk] . (3) 

Here e and rjj^k are random variables. The first one is 
a random number in the interval (0,1), which can be 
distributed either uniformly or with a certain probabil- 
ity distribution g{e), as in some generalizations of the 
basic Al-model [20]. The second one is a random di- 
chotomous variable responsible for the unidirectionality 
of the wealth flow as well as for the nonlinear charac- 
ter of the dynamics. It is a function of the difference 
between the wealths of the interacting agents j and k, 
rjj^k = (f'ixk — Xj), assuming the value ijj^k = 1 with 
probability po for Xj > Xfe or the value rjj^k = with 
probability 1 — po for Xk > xj . The value rjj^k = 1 pro- 
duces a wealth transfer ]Ax| = ecoxj from agent j to k, 
while the value rij,k = corresponds to a wealth transfer 
|Ax| = euxk from k to j. 
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(5) it follows that for n < 1, i.e., for A < 1/4 (w > 3/4), 
the F-distribution diverges for a; — > 0, as visible in Fig. 1 
for the cases A = 0, A = 0.1, and A = 0.2. For the critical 
value A — 1/4 {uj — 3/4), separating the distributions 
which diverge from those which go to zero for a: — > 0, an 
exponential distribution is obtained, 



(7) 



The Al-model has a simple mechanical analogue if 
the quantity D = 2n defined in Eq. (5) is interpreted 
as an effective dimension for the system and f3~^ as a 
temperature. It is easy to check that the distribution 
jn{l3x) = 7_D/2(/3a;) given by Eq. (4) is the equilibrium 
distribution for the kinetic energy of a perfect gas in D 
dimensions as well as for the potential energy of a D- 
dimensional harmonic oscillator or a general harmonic 
system with D degrees of freedom. This definition of 
effective dimension is consistent with the equipartition 
theorem, since 



(x) =nl3~^ = D/3-V2, 



(8) 



see Ref. [381 for details. 



B. A2-model 



FIG. 1: Equilibrium wealth distribution for the basic version 
of the Al-model defined by Eqs. (3), for the case po = 1/2: re- 
sults of numerical simulations (symbols) and fitting functions 
Eqs. (4) (curves) for different values of the saving parame- 
ter A = 1 — a; in linear (above) and semi-log (below) scale. 
The value of n is governed by Eq. (5). In this simulation the 
average wealth is (x) — 1. 



A special case of the Al-model is that with symmet- 
rical interaction, obtained for po — 1/2. Notice that for 
this value of po the random variable rjj^k = rj becomes 
independent of Xj and Xk ■ We have studied this particu- 
lar case through numerical simulations for various values 
of the saving parameter A. The system considered was 
made up of = 10^ agents with equal initial wealths 
— 1 and the transactions were performed until equi- 
librium was reached. The equilibrium distributions in 
Fig. 1 were obtained by averaging over 10^ different runs. 
They are well fitted by the F-distribution 



f{x) = Mn{Px) 



JL 



(/3a:)"~i exp(-/3x) , (4) 



where 



{x)/n, 
D _ 1 

y ^ 2(1 - A) 



2A 



3 

2uj 



(5) 
(6) 



Since A = l—uj £ [0, 1), the parameter n is a real number 
in the interval [1/2, oo). Notice that from Eqs. (4) and 



The One-Parameter Inequality Process model, here re- 
ferred to as A2, is another model introduced by John 
Angle and is described in detail in Refs. [22, 23]. It dif- 
fers from the Al-model considered above in that it only 
employs a stochastic dichotomic variable 'qjk, which can 
assume randomly the values rjjk = or rijk = 1. The 
model is defined by Eqs. (2) with 



Ax 



-rjjk^^Xk + (1 - rijk)i^Xj 



(9) 



The model describes a unidirectional flow of wealth from 
agent k toward agent j for 7]jk = 1 or vice versa for 
Vjk = 0. For the particular case in which the two val- 
ues of rjjk are always equiprobable, one can rewrite the 
process, without loss of generality, with a Ax = ujXk in 
Eqs. (2). Numerical simulations of this model confirm 
the findings of Refs. [22, 23], that for small enough oj the 
stationary wealth distribution is well fitted by a F-distri- 
bution 7„(a;), with n « 1/w — 1 = A/(l — A). We find 
that this fitting (not shown) is very good at least up to 
A « 0.7. 



C. B-model 

Another KWEM was introduced in 1988 by Eleonora 
Bennati [24, 25]. Its basic version, that we discuss here, 
is a simple unidirectional model where units exchange 
constant amounts Axq of wealth [24-26]. In principle, in 
the B-model a situation where the wealths of the agents 
would become negative could occur. This is prevented 
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FIG. 2: In the B-model the quantity x can only vary by a 
constant amount Aa; = Aa;o, which can e.g. be lost by a unit 
j and then absorbed by a unit k, analogously to the emission- 
absorption process of light quanta of constant frequency. 

allowing the transaction to take place only if the condi- 
tion x'j,x'i^ > is fulfilled, i.e., the process is described 
by Eq. (2) with Aa; = Axq if x^-, > and with Aa; = 
otherwise. Since the wealth can vary only by a con- 
stant amount Axq, the model reminds a set of parti- 
cles exchanging energy by emitting and re-absorbing light 
quanta, as illustrated symbolically in Fig. 2. Analytically 
the equilibrium state of the B-model is well described 
by the exponential distribution (7). A main difference 
respect to the other models considered here is that in 
the B-model the amount of wealth exchanged between 
the two agents is independent of Xi, while in the other 
models represents a multiplicative random process, since 
Ax oc Xi. 

D. C-model 




X 



FIG. 3; Equilibrium wealth distributions in linear and semi- 
log scale for different values of the saving parameter A in the 
closed economy model defined by Eqs. (10). The continuous 
curves are the fitting functions, i.e. a F-distribution of the 
form Eq. (4) just as for the Al-model, but the value of n(A) 
is given by Eq. (12). 



In the model introduced in 2000 by A. Chakraborti 
and B. Chakrabarti [27] the general exchange rule reads. 



Xxj + e (1 - A)(xj + Xk) , 
Xxk + e (1 - X){xj + Xk) , 



(10) 



where e = 1 — e. Here the new wealth x'^ (xj.) is expressed 
as a sum of the saved fraction Ax^ (AxJ.) of the initial 
wealth and a random fraction e (e) of the total remaining 
wealth, obtained summing the respective contributions of 
agents j and k. Equations (10) are equivalent to Eqs. (2), 
with 

Ax = Loilxj — exfc) = (1 — A)(exj — exfc) . (11) 

Like in the Al-model, at equilibrium the system is well 
described by a F-distribution (4). For the parameter n 
we find now [38, 39] 



D 1 + 2A 



n — — = 

2 1 - A 



(12) 



which is twice the value of the corresponding parameter 
of the Al-model with po — 1/2, discussed in Sec. Ill A. 



In Fig. 3 numerical results are compared with the fit- 
ting based on Eq. (12). In this case the probability den- 
sity is always finite for x — *■ 0, since for A = (a; = 1) one 
has n = 1 and the distribution does not diverge, being 
equal to the exponential function (7). 



E. D-model 

The models introduced in 2000 by A. Dragulescu and 
V. M. Yakovenko [28] were conceived to describe flow 
and distribution of money. They have a sound interpre- 
tation both of the conservation law x^ -)- xj. — Xj + Xk, 
since money is measured in the same unit and conserved 
during transactions, and of the stochasticity of the up- 
date rule, representing a randomly chosen realization of 
trade. Various models were considered in Ref. [28], with 
a Ax either constant (similarly to the B-model discussed 
above) or dependent on the values Xi of the agents; also 
more realistic models, in which e.g. firms were intro- 
duced or debts were allowed. For simplicity, we consider 
among them the model which probably best represents 
the random character of KWEMs, referred to as the D- 
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FIG. 4: The D-model as well as the C-model prescribe a mi- 
croscopic interaction between two units analogously to a ki- 
netic model of gas in which, during an elastic collision, two 
generic particles j and k exchange an energy amount Ax. 



model below, in which the total initial amount Xj + Xk is 
reshuffled randomly between the two interacting units, 



e {xj + Xk) , 
e {xj + Xk) ■ 



(13) 



Equivalcntly, the dynamics can be described by Eqs. (2), 
with 



Ax = 



exj - exk ■ 



(14) 



The D-model is formally recovered from the C-model for 
A = (w = 1). 

The equilibrium distribution of the D-model is well 
fitted by the exponential distribution (7). A mechan- 
ical analogue of the D-model is a gas, in which parti- 
cles undergo pair collisions in which some energy is ex- 
changed [40], as symbolically illustrated in Fig. 4. 



F. Stationary wealth distributions 

The parameters of the F-distribution, obtained from 
the fitting of the wealth distributions of the stationary 
solutions for the models considered, are summarized in 
Table F The analytical forms of the respective param- 
eters n, given as a function of uj or A, provide a good 
fitting: for the model A2, the fitting is good only up to 
A w 0.7. The close analogies among the various models 
are evident, however the existence of a general solution 
has not been demonstrated, see e.g. Refs. [23, 41]. 



TABLE I: Comparison of dependence of the fitting parameter 
n in the F-distribution (4) on A or oj, for the basic homoge- 
neous versions of the A-, C-, and D-models. 



Model 




) n(A) 




Al 


3/2uj- 


-1 (1-F2A)/2(1- 


-A) 


A2 


1/lo- 


1 A/(l-A) 




C 


3/lo- 


-2 (1 + 2A)/(1- 


A) 


D 


1 


1 





IV. INFLUENCE OF HETEROGENEITY 

Here we discuss the influence of heterogeneity, consid- 
ering as an example the generalization of the C-model. 
Heterogeneity is introduced by assigning a different pa- 
rameter uji (Xi) to each agent i. The formulation of the 
heterogeneous models can be straightforwardly obtained 
from those of the corresponding homogeneous ones by 
replacing the generic term uxi (Xxi) with UiXi (XiXi) in 
the evolution law. In the case of the C-modcl Eqs. (10) 
become 

x'j = XjXj + e [(1 - Xj)xj + (1 - Xk)xk] , 
x'k = XkXk + - Xj)xj + (1 ~ Xk)xk], (15) 
and the exchanged amount of wealth in Eqs. (2) is now 
Aa; = etOj Xj—euk Xk — €{1 — Xj)xj —e{l — Xk)xk ■ (16) 

The set of parameters {oJi} {{Xi}) is constant in time 
and specifies the proffles of the agents. The values {cui} 
{{Xi}) are assumed to be distributed in the interval be- 
tween and 1 with probabilities hi (gi) and ^ihi — 1 
(X^i 5i — !)■ the limit of an infinite number of agents, 
one can introduce a probability distribution h{Lo) [g[X)\, 
with jldujh{uj)^l [^ldXg{X) = \]. 

Various analytical and numerical studies of this model 
have been carried out [2, 19, 22, 36, 41-48] and as a 
main result it has been found that the exponential law 
remains limited to intermediate x-values, while a Pareto 
power law appears at larger values of x. Such a shape 
is prototypical for real wealth distributions. Numerical 
simulations and theoretical considerations suggest that 
the power law exponent is quite insensitive to the details 
of the system parameters, i.e., to the distribution h{uj). 
In fact, the Pareto exponent depends on the limit g{X^ 
1). If g{X) ~ (1 - A)""i with A ^ 1 and a < 1, then the 
corresponding power law has an exponent a [44]. Thus, 
in general, agents with A^ close to 1 are responsible for 
the appearance of the power law tail [44, 47, 49]. 

Probably the most interesting feature of the equilib- 
rium state is that while the shape of the wealth dis- 
tribution fi{x) of agent i is a F-distribution, the sum 
of the wealth distributions of the single agents, f{x) = 
'^ifi{x), produces a power law tail. Vice versa, one 
could say that the global wealth distribution f{x) can be 
resolved as a mixture of partial wealth probability den- 
sities fi{x) with exponential tail, with different param- 
eters. For instance, the corresponding average wealth 
depends on the saving parameter as (a;)i oc 1/(1 — AJ = 
l/uji] see Refs. [47, 49, 50] for details. 

Importantly, all real distributions have a finite cut- 
off; no real wealth distribution has an infinitely extended 
power law tail. The Pareto law is always observed be- 
tween a minimum wealth value Xmin and a cutoff a;inaxj 
representing the wealth of the richest agent. This can 
be well reproduced by the heterogeneous model using an 
upper cutoff Amax < 1 for the saving parameter distribu- 
tion g{X): the closer to one is Amax, the larger is Xmax and 
wider the interval in which the power law is observed [49] . 
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The role of the A-cutoff is closely related to and rel- 
evant for understanding the relaxation process. The re- 
laxation time scales of single agents in a heterogeneous 
model are proportional to 1/(1— Aj) [37]. This means that 
the slowest convergence rate is determined by 1 — Amax- 
In numerical simulations of heterogeneous KWEMs, one 
necessarily employes a finite A-cutoff. However, this 
should not be regarded as a limit of numerical simula- 
tions but a feature suited to describe real wealth dis- 
tributions. Simulations confirm the fast convergence to 
equilibrium for each agent with the above mentioned time 
scale [37]. Gupta has demonstrated numerically that the 
convergence is exponentially fast [51]. 

In Ref. [52] it has been claimed that heterogeneous 
KWEMs with randomly distributed Aj (0 < Ai < 1) 
cannot undergo a fast relaxation toward an equilibrium 
wealth distribution, but the relaxation should instead 
take place on algebraic time scale. This in turn means 
that there cannot exists any power law tail. Such claims 
are probably correct for systems with a A-distribution 
g{X) rigorously extending as far as A = 1, corresponding 
to a power law tail extending as far as a; = oo. However, 
this does not apply to KWEMs with a saving parameter 
cutoff Amax < 1, which is the natural choice in describing 
real systems, as well as in numerical simulations, employ- 
ing a finite A-cutoff: in this case the largest time scale is 
finite and relaxation is fast. 



V. GENERALIZATIONS 

In this section a unified reformulation of the exchange 
laws of the A-, C-, and D-models is suggested and as an 
example an application to the C-model is made. 

A. Reformulation 

It is possible to reformulate the evolution law either 

through a single stochastic saving variable A or an equiv- 
alent stochastic exchange variable uj ~ I — X. This for- 
mal rearrangement of the equations maintains the form 
of the evolution law very simple and has at the same time 
the advantage to be particularly suitable to make further 
generalizations. For the sake of generality, wc consider 
the case of a heterogeneous system characterized by a 
parameter set {u)i}. The models discussed above (apart 
from the B-model) can be rewritten according to the ba- 
sic equations (2), where the wealth exchange term is now 
given by 

Aa; = Uj Xj — UkXk ■ (17) 

The meaning of the new stochastic variables Qj and u>k 
introduced is simple: ujj represents the fraction of wealth 
given by agent j to k during the transaction, and vice 
versa for Wk- Comparison with the equations defining the 
A-, C-, and D-models provides the following definitions 
for u)j and Qk- 



• In the Al-model, ujj and u)k are independent non- 
linear stochastic functions of the agent wealths Xj 
and Xk, 

i^k = {'i--Vj,k)eu)k, (18) 

where r/j^k = <f>{xk — Xj) = 1 with probability po for 

Xk — Xj > and rjj^k = with probability 1 — po for 
Xk — Xj < 0, while e is a random number in (0, 1). 
For rjj^k = one has Qj = and ujk G (0,0;^), 
whereas for r]j k = 1 one has Uk = and Uj € 
(0,a;,). 

• In the A2-model, Uj and Uk only contain the di- 
chotomic variable, 

u)j = rij,k'^j , 

Wfe = (1 - Vj,k)i^k ■ (19) 

• For the C-model, 

ojj = eojj , tJj e (0, CiJj) , 

Wfe = (l-e)cL)fe, (I;fee(0,Wfe), (20) 

where e is a random number in (0, 1). 

• The D-model is recovered from Eqs. (20) of the C- 
model when u)i = 1 for each agent i. 



TABLE II: Comparison of the A-, C-, and D-models: explicit 
forms of the random variables Cjj and Qk in the unified refor- 
mulation (17). See text for details. 



Model Qj uJk 



Al 




{l-rj)eu}k 


A2 




{l-rj)iJk 


C 




(1-r) '^k 


D 


(' 


Ci-r) 



The reformulation is summarized in Table II with ref- 
erence to Eq. (17). Different generalizations can now be 

done changing only the properties of the stochastic vari- 
ables LOi, while maintaining the same formulation (17) of 
the exchange law. 

B. An example 

As an example which can be represented through 
Eq. (17), we consider a generalization of the homoge- 
neous C-model. In the original version there is a con- 
straint on the maximum fraction of invested wealth, given 
by a value < a; < 1 of the exchange parameter, or 
equivalently on the minimum saved fraction, given by a 
value < A < 1 of the saving parameter. Now an addi- 
tional constraint on the minimum fraction of the invested 
wealth is assumed. This may describe e.g. trades which 
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FIG. 5: Fitting parameter n(A, A') of tiie F-distribution 771(2;) 
for the generalized model with A < Ai < A'. Dotted and 
dashed lines: n as a function of A (< A' by definition) for 
fixed values of A' shown on the right side. These lines stop at 
A = A' since by definition A < A'. Continuous lines: n as a 
function of A' ((> A) for the same fixed values of A shown on 
the right side. 



always have a minimum risk for an agent. It can be rep- 
resented by an analogous parameter lu' , with < uj' < to, 
representing the minimum fraction of wealth invested in a 
single trade. One can also define a parameter A' — l~uj', 
with A < A' < 1, representing the maximum fraction of 
saved wealth (i.e. it is not possible to go through a trade 
without risking a non-zero amount of wealth). Then the 
stochastic variables Ui in Eq. (17) become uniform ran- 
dom numbers in intervals defined by the parameters lu' 
and Lo (or by A' and A), 



G(c.',c^) = (l-A',l-A), 



(21) 



We have performed numerical simulations for a set of 
combinations of parameters (A, A') and found that the 
equilibrium distributions are always well fitted by the 
same T-distribution (5). However, we have not found a 
simple analytical formula for fitting the dependence of 
the parameter n on the saving parameters A and A'. The 
behavior of n versus A (A') is represented graphically in 
Fig. 5. Dotted/dashed curves (different colors) represent 
n versus A for the different fixed values of A' shown on the 
right side These curves stop at A = A', since by definition 
A < A'. From there the continuous (red) curves start, 
which represent n versus A' for the same fixed values 
of A listed in the legend on the right. The first (dashed 
green) curve from the top extending on the whole interval 
A = (0, 1) represents n as a function of A f or A' = 1 
and corresponds to the original homogeneous C-model. 
For this particular case n{X) is known to diverge as n ~ 
1/(1 - A) - l/w for A ^ 1 (see Table I), while in all the 
other cases n is finite. 



VI. CONCLUSIONS AND DISCUSSION 

We have reviewed some basic KWEMs of closed econ- 
omy systems, introduced by scientists working in differ- 
ent fields, allowing us to point out analogies and differ- 
ences between them. We have first considered the ho- 
mogeneous models and then discussed the influence of 
heterogeneity. The heterogeneous KWEMs are particu- 
larly relevant in the study of real wealth distributions, 
since they can reproduce both the exponential shape at 
intermediate values of wealth as well as the power law 
tail. 

In all the models discussed, including the heteroge- 
neous one, the equilibrium wealth distribution of a single 
agent is well fitted by a F-distribution, known to be the 
canonical distribution of a general harmonic system with 
a suitable number of degrees of freedom. This suggests a 
simple mechanism underlying the (approach to) equilib- 
rium of these systems, similar to the energy redistribution 
in a mechanical system. However, a general demonstra- 
tion that the F-distribution is the stationary solution of 
KWEMs and an understanding of how it arises is still 
missing (see Refs. [23, 41, 53, 54] for theoretical consid- 
erations on and the microcconomic formulation of this 
issue). 

Furthermore, we have discussed how in a heteroge- 
neous KWEM the sum of the single agent wealth dis- 
tributions can produce a power law tail. In particular, 
we have clarified some issues concerning the relaxation 
process and the existence of power law tails: whenever 
there is a finite cutoff in the saving parameter distribu- 
tion, the largest time scale of the system is finite and one 
observes a fast (exponential) relaxation toward a power 
law, which extends over a finite interval of wealth. The 
width of such interval depends on saving parameter cut- 
off. 

Due to the similarity of the structures of the models 
discussed, we have proposed a novel unified reformulation 
based on the introduction of suitable stochastic variables 
uji, representing the actual fraction of wealth lost by the i- 
th agent during a single transaction. This unified formu- 
lation lends itself easily to further generalizations, which 
can be obtained by modifying the stochastic properties of 
the variable LUi only, while leaving the general evolution 
law unchanged. We have illustrated the new formulation 
by working out in detail an example, in which the frac- 
tion of wealth lost Coi is characterized by a lower as well 
as an upper limit. 

Besides the KWEMs considered in the present pa- 
per, originally formulated through finite time differ- 
ence stochastic equations, other relevant (versions of) 
KWEMs have been introduced in the literature; see 
Refs. [10, 19, 55] for an overview. Their mathemati- 
cal formulation can be similar to the one of the present 
paper [21-23, 36, 56], or different approaches can be 
used, such as matrix theory [57], the master equa- 
tion [35, 58, 59], the Boltzmann equation [41, 52, 60-63], 
the Lotka-Volterra equation [64, 65], or Markov chains 
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models [66 68] . All these models share a description of 
wealth flow as due to exchanges between basic units. In 
this respect, they are all very different from the class 
of models formulated in terms of a Langevin equation 
for a single wealth variable subjected to multiplicative 
noise [29, 69-72]. The latter models can lead to wealth 
distributions with a power law tail. In fact, they converge 
toward a log-normal distribution, which, however, does 
not flt real wealth distributions as well as a F-distribution 
or a /3-distribution and is asymptotically characterized by 
too large variances [21]. 

Finally, we would like to point out that even though 
KWEMs have been the subject of intensive investiga- 
tions, their economical interpretation is still an open 
problem. It is important to keep in mind that in the 
framework of a KWEM the agents should not be related 
to the rational agents of neoclassical economics; an in- 
teraction between two agents does not represent the ef- 
fect of decisions taken by two economic agents who have 
full information about the market and behave rationally 
in order to maximize their utility. The description of 
wealth flow provided by KWEMs takes into account the 
stochastic element, which docs not respond by defini- 
tion to any rational criterion. Also some terms employed 
in the study of KWEMs, such as saving propensity (re- 



placed here by saving parameter), risk aversion, etc., can 
be misleading since they seem to imply a decisional as- 
pect behind the behavior of agents. Trying to interpret 
the dynamics of KWEMs through concepts taken from 
the neoclassical theory leads to obvious misunderstand- 
ings [73]. However, it is interesting to note that very 
recently, Chakrabarti and Chakrabarti have put forward 
a microeconomic formulation of the above models, using 
the utility function as a guide to the behavior of agents in 
the economy [54]. Instead, KWEMs provide a descrip- 
tion at a coarse grained level, as in the case of many 
statistical mechanical models, where the connection with 
the microscopic mechanisms is not visible; however, the 
equivalence is maintained. 
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